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Abstract 



' In this article we consider the system of equations Au; = pi (x) fi (ui, for i — 1, on 

, N > 3 and d £ {1,2,3,4, ...}. We prove that the considered system has a bounded positive 
entire solution under some conditions on pi and fi. Also, we give a necessary condition as well as 
a sufficient condition for a positive radial solution to be large. The method of proving theorems is 
essentially based on a successive approximation. Furthermore, a non-radially symmetric solution is 
obtained by using a lower and upper solution method. 
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O ■ 1 Introduction 



In this paper we study the existence of solutions for the semilinear elliptic system 

Ami = pi (x) fi (ni, ...,Ud) , 

inM^, (1.1) 

Aud = Pdix) fd{ui,...,Ud) , 

where A'' > 3, d > 1, and the functions pj, fj (j = are supposed to satisfy the following 

hypothesis: 

(PI) Pj : — 7- [0, oo) are locally Holder continuous functions of exponent a G (0, 1); 

(CI) fj : [0,00)'^ — )• [0,00) are continuously differentiable in each variable, /j(0, ..,0) = 0, and 

fj (si, Sd) > if Sj > for some i = 1, . . . , d; 

(C2) fj are increasing in each variable; 

(C3) /,°°[F(.)]-V2d. = oo {F{s) = f^^§J,{t,...,t)dt). 

Similar problems to those we are analyzing here are related to steady-state reaction-diffusion, 
subsonic fluid flows, electric potentials of some bodies and control theory. 
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As an example, the first motivation for studying the above problem stems from the article [H] where 
the reader observe that such problems arise from the description of the basic stochastic control theory. 
The controls are to be designed so that the state of the system is constrained to some region. Finding 
optimal controls is then shown to be equivalent to finding large solutions for a second order semilinear 
elliptic partial differential equation. In terms of the dynamic programming approach, an explosive 
solution of (jl.ip corresponds to a value function (or Bellman function) associated to an infinite exit 
cost (see [8]). 

Another motivation comes from the work of [3] where the parabolic problem corresponding to 
system (jl.ip are models of steady state of non-linear heat conduction through a 2-components mixture. 

The numerous applications that lead to favorable answers in order to establish new and significant 
results for problem (jl.ip as well as the recent results in the field motivate the study of more generally 
class of problems ()l.ip . 

The main results of this paper are the following two theorems. They substantially solve the open 
problem proposed in [7j and complete the results of [2] where only sufficient conditions are obtained. 

Theorem 1.1. Suppose that (PI), (C1)-(C3) are satisfied. If there exists a positive number e such 
that 

t S (p,; it)dt < oo where 09„- (t) = maxp,(a;), (1-2) 



d 

and r S (r) is nondecreasing for large r, then system has a nonnegative nontrivial bounded 
j=i 

solution on M^. //, on the other hand, pj satisfy 

d 

I t T, 1/;^ (t) dt = oo where t/j^ (t) = minp, (x) , (1-3) 
Jo J=i kl=* 



and r^^ ^ 
bounded radial solution on M^. 



S (r) is nondecreasing for large r, then system has no nonnegative nontrivial entire 



Theorem 1.2. Assume that pj : M.^ — )■ [0,oo) (j = l,...,d) are spherically symmetric continuous 

functions {i.e. pj (x) = pj {\x\)). If fj {j = 1, . . . ,d) satisfy (CI) -(C3), then the problem has a 

d 

nonnegative nontrivial entire radial solution. Suppose furthermore that r^^~^ S pj (r) is nondecreasing 
for large r. If Pj satisfies 



t 

s^~^Pj (s) dsdt = oo for all j = 1, . . . , d, (1-4) 



then any nonnegative nontrivial solution {ui,...,Uct) of is large. Conversely, if U.l\) has a non- 

negative entire large solution, then pj satisfy 



oo d 

r^+^ S Pj (r) dr = oo, (1.5) 



for every e > 0. 
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2 Preliminary result 



We recall the following definition of lower and upper solution which are our main tool in the proof of 
the solvability of problem (II. ID . 

d 



C, 



2,a 
loc 



Definition 1. A function [wi, ...,Wd) G 
problem if 

Awi>Pi{x) fi {wi, ...,Wd) in M.^ for all i = l,...,d. 



(a G (0, 1) ) is called a lower solution of the 



Definition 2. We say that {vi, ■■■,Vd) G C^j" (M ) (a G (0, 1)) is an upper solution of the problem 

onp if 

Avi<pi{x) fi{vi,...,Vd) inR^ for all i = 1, ...,d. 

We need the following lemma which can be found in [Jl Theorem 5.1, pp. 146]: 

Lemma 2.1. Make the same assumptions on pj and fj (j = l,...,d) as in Theorem li.il If the 
problem li.i)) has a pair of upper and lower bounded solutions {vi,...,Vd) and {wi, ...,Wd) fulfilling 
Wi{x) < i = l,...,d,\/x G then there exists a bounded function {ui,...,Ud) belonging to 



CI 



2,a 



(a G (0, 1) ) with 

Wi{x) < Ui{x) < Vi{x), i = 1, d, Vx G 



and satisfying 



3 Proofs of main results 

In this section we give the proofs of Theorems 11.11 and 11.21 



Proof of the Theorem 11.11 Assume that ()1.2p holds. In this section, we use the method of 
upper and lower solutions to show the existence of positive bounded solutions of the problem (II. ip . 
Our aim is to construct an upper solution (ui, ...,Vd) and a lower solution {wi, .■.,Wd) for the problem 
(|l.ip . On this purpose, we first prove the existence of {wi^ ...,Wd) to the system 

Awx (r) = (j) fi (wi, ...,Wd) for r := |x| , 
i ... inR^. (3.1) 

Awd{r) = (Pd{r) fd{wi,...,Wd) for r := |x| , 

Observe that we can rewrite (|3.ip as follows: 

{r^~^wi (r)) = Vi (r) fi {wi, ...,Wd) , 
{r^"^Wd [r])' = r^-'^ip^ (r) fd (wi, .-^Wd) ■ 
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Then radial solutions of are any solution (wi, Wd) of the integral equations 



Wd{r) = 2 + lo W=T Jo Vd (s) /d (w'l (s) , ■■■,Wd (s)) dsdt. 



To establish a solution to this system, we use successive approximation. Define sequences <,w^} 

I ■>)j=i,...,d 

on [0, oo) by 



= i, r > 0, 



^ ir) = l + Jo W-T /o V. (^) (^r' (^) , -.w',-' is)) dsdt, i = 1, d. 
We remark that, for all r > 0, j = 1, d and k & N 

Moreover, proceeding by induction we conclude t^'? f are non-decreasing sequence on [0, oo). 

L -'J j=l,...,d 

r k^^-^ 

We note that iw^} satisfy 

I ) j=l,...,d 



By the monotonicity of \ Wj > we have the inequalities 

I J 7=1,. 



fc>i 



N-l f„„k\' 



""^"Vd (0 fd (wt\ ^'rf-') < r^-Vd (r) S /i f S (r) , S (r) 

\ / 1=1 \ 7 = 1 7 = 1 



Thus, summing up inequalities, we have 



1=1 \ 



< r^-\S (r) S /, ( S (r) , S (r) 
1=1 1=1 \j=i j=i 



(3.2) 



(3.3) 



Choose i? > so that r S • (r) are non-decreasing for r > R. We are now ready to show that 

.7 = 1 



Wj (R) and iwj (R)] , both of which are nonnegative, are bounded above independent of k. To do 



this, let 



max{ipj (r) : < r < R}, j = 1, d. 
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Using this and the fact that iw'-] > 0, we note that (13. 2p yields 



w 



<0f s/, ( ^wf,...,j:wl 



1=1 \i=i 



' i=i \i=i 



1=1 



<cp!lj:fi(j:wl...,j:w. 



1=1 



This imphes that 



d 

E 
1=1 



d ^\ d / d , d , 

, 1=1 I 1=1 \ 1=1 1=1 



On the other hand, multiplying this equation by ( ^^wf ) and integrating we also have 



1=1 



< 2 ( ^^(j^ij J .S /i (s, s)ds,0<r< R. 



(3.4) 



Integrating the above inequality, we see that 



i—l 



S fi{s,...,s) ds 



i=l 



dt < J2 ( \R. 

d 

It follows from the above relation and by the assumption C3) that S wf (R) is bounded above inde- 

i=l 

pendent of k. Using this fact in (j3.4p shows that the same is true of f S [R)] . Thus, the sequences 



i=l 



W. 



f (R) and (^wf (R))' are bounded above independent of k. Finally, we show that the non-decreasing 



sequences is bounded for all r > and all k. Multiplying (13. 3p by r y ) ^-^id integrating 

gives 



< 



r^-i(.Su;f(r) 



( (R) 



r d d / d , d \ / d \' 

+ 2 / 52A^-2 S y,, is) E h E wl S S ds, 

Jr «=i *=i V*=i / V=^ J 



for r > R. Noting that, by the monotonicity of s 5] 99^ (s) for s > R, we get 



( Ewf (r) 



< C + 2r2^-2 S (r)F(s mf(r 

4 = 1 \ '1 = 1 



where C 



i=l 



, which yields 



F ( (r) 



1/2 



(3.5) 
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and hence 



dr 



Integrating this and using the fact that 



d \ V2 



for each e > 0, we have 

d 



F ( ^w^ (r) 



1=1 



-1/2 



-1-e 



+ 2S (^, (r) 

2=1 



1/2 



1/2 



< S 



l+e 



-1-e 



.S«;f(r) 



r ,1-N 



i=l 



-1/2 



dt 



+ 



1/2 



< </C 



F ( .E (R) 



r^e-^dt 



+ 



The above relation is needed in proving the bounded of the function s > in the fohowing. 

Indeed, since for each e > the right side of this inequahty is bounded independent of k (note that 



(3.6) 



wf {t) > l/d for all i = 1, d), so is the left side and hence, in light of C3), the sequence < S w: 

i=l 



k>l 



is a bounded sequence and so <Wj > are bounded sequence. Thus tt^? f — > {wj}j=i a as 

I J 3=1, ■■■,d J j=l,...,d ■' 

A; — 7- oo and the limit functions {wj} -^^ ^ are positive entire solutions of system (jS.ip . It is obvious 
that every solution of ()3.ip is a lower solution to (jl.ip . 



We let M be the least upper bound of S Wi and note that 

i=l 



M = lim S (r) . 

r— >oo i=l 

Now let tp^ (t) = mini^i^iPi (x) {i = 1, d) and {t^jj^-^]^ ^ be the positive increasing bounded solutions 
of 

\,{r) = M + /; s^- Vi (s) /i (^;i (s) , va {s)) dsdt, 



Vd{r) = M + /q /q (s) /d (^'1 (s) , Vd (s)) dsfit, 

which, of course, satisfies (j3.ip with (j = 1, d) replaced with vj {j = 1, ...,d) and ifj (j = 1, d) 
replaced with ipj (j = 1, ...,d). It is also clear that > M. If we modify the proof of the lower solution 
{wi, Wd) we obtain the existence of an upper solution {vi, Hence the standard upper-lower 

solution principle implies that (jl.ip has a solution (ui, ■■■^Ud) such that 



i(^i(x) < Ui{x) < Vi{x), i = 1, d, 
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(see Lemma 12. ip . 

Now assume that ()1.3p holds. Assume to the contrary that there exist a nonnegative nontrivial 
entire bounded radial solution {ui, ■■■^Ud) on for the system (jl.ip . Assuming Mj = sup^gjgjv Ui (x) 
{i = 1, ...,d) and knowing that u'- > 0, we get liuir^ooUi = Mi. Thus there exists i? > such that 
Ui > Mi/ 2 (i = 1, ...,d) for r > R. From conditions of /j, it follows that 

h {ui, ...,Ud) > h {Ml/2, Md/2) := 4, for r > i? and i = 1, ...,d. 

Using this we get the following 

ui > ui{0)+ r t^-"" f s^-^iPi{s)fi{ui{s),...,Ud{s))dsdt 
Jo Jo 



(3.7) 



t 

l-N I „N-1 



Ud > Ud (0) + I t^-^ [ 
Jo Jo 



i^d (s) fd (ui (s) , Ud (s)) dsdt. 



Rearranging the terms, and by using these conditions in ()3.7p follows 

d I r d 

S Ui (r) > dull + m2 — / t T, ipi (t) dt ^ co as r ^ oo, 

i=l N Jj^ i=l 

where mi = minjui (0) , (0)} and 1712 := minjcQ, ...,Cq}. A contradiction to the boundedness of 

d 

S Ui (r). This proves Theorem ll.il 



1=1 



Proof of the Theorem 11.21 It is known from [6l [121 Theorem 2] that the problem 



Az (r) = Pi (r) S /j {z (r) , z (r)) for r := \x\ , x £ 

1=1 i=l 



has a non-negative non-trivial entire solution. Moreover, for each R > 0, there exists cr > such that 
z (R) < cr. Due to the fact that z is radial, we have 

r 1 AT ^ d d 

z{r) = z (0) + / -jT— / S Pi (s) S fi {z {s) , z (s)) dsdt for all r > 0. 

Jo t Jo 

We choose f3i G (0, z (0)]. Define the sequences \u^ r on [0, 00) by 

= ... = n[J = /3i for ah r > 

ui (r) = /?! + /; ^ Jo* V (s) /i {ut' is) , is)) dsdt, 

, n^(r) = /3i + /;pv^/o*s^-ip,(5)/rf(n^i(s),...,n^i(.))dsdt, d > 1. 
With the same arguments as in the proof of Theorem 1 1.1 1 we obtain that \u^ f are non-decreasing 

L J j=l,...,d 

sequence on [0, 00). Because z' (r) > follows < I3i < z (0) < z (r) for all r > and so 



1 



ulir) = Pi+ I I s^-'piis)f^{uiis),...,u^dis))dsdt 



< z{0) + 



fN-l , 
JO 

1 



Jo 



.N-1 



S Pi (s) S /j (z (s) , z (s)) dsdt = z (r) 



i=i i=i 
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Thus u\ (r) < z (r). Similar arguments show that 



Uj<z (r) {j = 1, d) for all r G [0, oo) and k > 1. 

Thus, for every x G the sequence Uj {\x\) {j = 1, ...,d) have subsequences, denoted again by 
(j = l,...,d), converging and 

(ni (|x|) , ...,Ud (|x|)) := ( lim u'l (|x|) , lim {\x\) 

\ fe— >00 K— >00 

is an entire radial solution of system (jl.ip . 

Let (ui, ■■■,Ud) be any non- negative non-trivial entire radial solution of (jl.ip and suppose that pj 
(j = 1, ...,d) satisfies (II. 4p . Since Ui [i = 1, is nontrivial and non-negative, there exists i? > so 
that Ui (R) > 0. Since u'- > 0, we get Uj (r) > Uj (i?) for r > R and thus from 

ui (r) = (0) + /q V (s) fi {ui (s) , Ud {s)) dsdt, 

Ud {r) = Ud (0) + /q s^~'^pd (s) /d (ui (s) , Urf (s)) dsdt, 

we get 



(r) = Ui (0) + /J' /g (s) /i (ni (s) , -u^ (s)) dsdt, i = 1, d 

> Ui {R) + fi (ui (i?) , -Urf (R)) s^~-^Pi (s) (isdt cx) as r oo. 

Conversely, if fi {i = l,...,d) satisfy (C1)-(C3) and {wi,..,Wd) is a nonnegative entire large solution 
of (jl.ip . then Wi satisfy 

[r^^i =pi (x)/i (u;i,...,u;d) , 

[r^^i (u;^)']' =Pd(2;) fd{wi, ...,Wd) ■ 
d 

Then, using the monotonicity of S pj (r) we can apply similar arguments used in obtaining 

Theorem 1 1.1 1 to get 



{^iw^ (r)) < VCr^~^ + J 2i^<f, (r) J F (^iw^ (r)) 



This conclusion and relation ()3.6p imply that 



Z — 1 



<CR + /^ti+^.Spi(t)dt, 



(3.9) 
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where 



Cr 



'C 



F\Y.w, {R) 



-1/2 ^Ar_2 



N 



+ 



By taking r — t- oo in (j3.9p we obtain that pj (j = 1 , . . . , d) satisfies (jl.Sp . This completes the proof of 
the theorem. 

We conclude this paper with some remarks: 

Remark 1. The time-independent Schrddinger equation in quantum mechanics is 

(/iV2m)An = {V - E)u 

where h = 6.625 • 10~^^ erg sec is the Planck constant, m is the mass of a particle moving under the 
action of a force field described by the potential V{x, y, z) whose wave function is u{x, y, z, t) and the 
quantity E is the total energy of the particle, problems which falls into the class of equations discussed 
here. 



Remark 2. // (C1)-(C3) are satisfied then 

fi (s, ••■,s) 
Remark 3. (see U^) If Cl)-C2) and 



-1/2 



dsdt = oo for all i = 1, .., d. 



-1 



oo / 

S fi ) ds = oo, 



rco / r-t d \ -1/2 

j ij .S^/i (s, ...,s)J dsdt = oo. 



are satisfied, then 
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